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We present a new approach to real-space multiple-scattering theory for molecules and clusters, 
based on the two-potential (distorted-wave) Lippmann-Schwinger equation formalism. Our approach 
uses a recently developed form [D. L. Foulis, Phys. Rev. A 70, 022706 (2004)], for the partial-wave 
expansions of the exact time-independent single-particle Green function for a general potential, to 
solve exactly the scattering problem for the distorting potential. The multiple-scattering problem 
for the full multicenter molecular potential is then developed along familiar lines, within a partition 
of space consisting of non-overlapping atomic spheres, but relative to the distorting potential. To 
achieve this some new general Green-function re-expansion formulas are derived, as well as further 
developments of our earlier partial-wave expansions. In passing, we make use of some of our results 
to derive a complete solution to the single-scattering problem for an arbitrary noncentral potential, 
and investigate some of the algebraic properties of Gaunt coefficients which arise frequently in 
our formulas. Based on the division of the multicenter molecular potential into the non-singular 
distorting potential and a remaining singular part we develop explicitly the secular equations of 
our approach and prove a result concerning the symmetry of the atomic matrices. The new secular 
equations are similar in overall form to those of related methods, but do not require any volume 
integrals. In our approach the computational burden consists essentially of the solution of the 
coupled radial Schrodinger equations for each atomic center, once in the actual atomic potential, 
and twice (giving the regular and irregular solutions) in the distorting potential; followed by integrals 
of combinations of these solutions over the atomic-sphere surfaces. We comment on key aspects of 
the numerical implementation, notably issues related to the choice of distorting potential. We treat 
both continuum (scattering) states and bound states within the same framework, and consider also 
the case of an outer sphere. 

PACS numbers: 03.65.Nk 31.15.A- 34.80.Bm 



I. INTRODUCTION 

In previous work [l[ the author has derived partial- wave expansions for the exact time-independent Green function 
of a general noncentral one-electron potential in terms of matrix solutions of the coupled radial Schrodinger equations 
for that potential. Although hopefully of interest in its own right to quantum-scattering theorists, the motivation 
was, in fact, to provide a key result in a new approach to multiple-scattering theory for general potentials, which is 
intended to be less cumbersome, more transparent, and more amenable to accurate numerical implementation than 
alternative approaches. In this article the full development of our new approach is presented. 

We consider here multiple-scattering theory in the restricted sense of the well known multiple-scattered- wave (MSW) 
theory of Johnson, Slater and co-workers (see for example Ref. fH) for bound states, extended by others [|| to 
continuum states and scattering, as well as photoionization [4]. It is expected, nevertheless, that some of the results 
developed below will have wider application. The MSW theory is itself an adaptation to finite clusters and molecules 
of the famous KKR method from solid-state physics. Although we concentrate here on methods for finite systems, 
where the practical applications are varied and numerous, it is expected that the approach of this present work may 
be adapted to the original case of periodic systems. 

As is well known, the main limitation for the original MSW and KKR methods was the requirement of the so called 
"muffin-tin" (MT) approximation for the potential. In this approximation space is partitioned into disjoint atomic 
spheres and a remaining interstitial region. The atomic potentials are taken to be spherically symmetric and the 
interstitial potential constant, so that the motion of an electron in the potential may be described in terms of free 
propagation between distinct atomic scattering events, for which the associated T matrices are diagonal. There have 
been considerable efforts over many years to extend these multiple-scattering methods to more general potentials, 
which began to bear fruit towards the end of the 1980s, and have led to some computational implementations. 

One of the main reasons that so much work has been done in this area, despite the fact that a variety of alternative, 
non-scattering methods have become available (in certain fields of application) , is that few, if any, of these alternatives 
have all the attractive features of the MSW/KKR approach. In particular, these are real-space, non-basis-set methods 
which combine a rapidly convergent partial-wave representation of multicenter wavefunctions with an accurate treat- 
ment of atomic cores. Furthermore, both bound and continuum states may be treated within the same framework, 
making the methods applicable to a wide range of problems. An attractive feature (of the original formalism at 
least), both from a conceptual and a practical point of view, is the separation between geometry and scattering that 
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is manifest in the secular equations of the method. 

The approach we shall present here has roots in the work of Natoli, Benfatto, and Doniach [5| who developed 
a particularly clear formulation of MSW theory for general potentials, applied to molecules and clusters, drawing 
themselves on the work of Beleznay and Lawrence Q, and Lloyd and Smith Q. In their approach the atomic T 
matrices are generalized to describe scattering from a non-spherical atomic potential, achieved by integrating the 
coupled radial Schrodinger equations in matrix form. The resulting T matrices then have non-zero off-diagonal 
elements. The effect of a varying interstitial potential on the electron scattering is taken into account via an extra 
term in the MSW secular matrix, related to a T matrix for the interstitial potential. 

This full-potential (FP) MSW method was implemented by the present author and colleagues [1] as the suite 
of computer codes fpx, and tested for analytically known bound-state and continuum problems as well as realistic 
molecular systems [t| [in the context of ab initio calculations in X-ray absorption fine-structure spectroscopy (XAFS)]. 
The codes have proved a useful pathfinder in exploring the benefits of full-potential calculations and have produced 
results [Io| which remain a benchmark pH in the field. They were designed to accommodate arbitrary molecular 
geometries, although optimized for symmetric systems, implementing the algebra of Ref. 5] completely, except for 
the interstitial T matrix which is calculated in the Born approximation. They have also had a useful general role in 
clarifying the relative importance of elements in the hierarchy of approximations underlying XAFS calculations. 

In practice, the fpx codes have proved of limited applicability, primarily because of the Born approximation for the 
interstitial T matrix, and the associated scaling behavior resulting from the volume integrals. A scheme to calculate 
the full interstitial T matrix was developed by Wastberg [13], and implemented for bound states of diatomic molecules. 
It is also a volume integral method relying on the use of an accurate quadrature grid for the interstitial region. Had 
development of the fpx codes been continued his scheme would undoubtedly have been incorporated. 

Other schemes to extend MSW/KKR methods to non-MT potentials have been developed, notably the "cellular 
method" in the context of band structure calculations for periodic systems, proposed by Williams and Morgan [l3j |. 
developed in a variety of forms and worked on by many others. (To find a way into the literature of this approach see 
Refs. Ul QUI] an d further references therein.) 

Despite a measure of progress, represented by the above approaches (among others), it appears still to be the case 
that non-MT MSW/KKR methods have been less widely taken up than might be expected. It is true, in general, 
that all of these extensions are significantly more complicated algebraically than the original MT-based methods, and 
their numerical implementation correspondingly more challenging. We shall not consider here in detail the reasons 
for the poor take-up of these methods, although it seems likely to this author that they are not unrelated to the latter 
observations. 

It is perhaps not immediately obvious how to improve the situation. The experience of the present author suggests 
that the emphasis should be on reducing the computational burden and, if not reducing algebraic complexity, at least 
increasing the transparency of the theoretical approach. Also crucially important are the details of the numerical 
implementation; although these will, of course, be dictated to a great extent by the theory. This present work is 
an attempt to develop a theory with reduced computational requirements, which is also well adapted to an accurate 
numerical implementation. 

Our particular approach has grown from the simple notion, mentioned in Ref. [l|, that a Green function better 
adapted than the free Green function, to the case of an electron propagating in a varying potential, should play 
a central role in the theory. This leads naturally to a consideration of the distorted-wave Lippmann-Schwinger 
equation and how it may be adapted to the multiple-scattering problem, and thus to the development of the necessary 
mathematical machinery. As will be seen below the resulting method retains the partition of space with distinct 
spherical atomic volumes, but requires no volume integrals. The computational burden is mainly the integration of 
coupled radial Schrodinger equations for each atomic center, and the evaluation of surface integrals over the atomic 
spheres. It is hoped that the algebraic development may be considered transparent and convincing, and perhaps 
conceptually simpler than other approaches in this field. 

II. THEORY 

We consider a molecule or finite cluster of N atoms with nuclear positions given by the vectors (for i = 1, . . . , N) 
relative to some origin. Associated with the molecule is a one-electron potential V(r), arising perhaps in the context of 
Hartree-Fock or density-functional calculations, of a multicenter character, i.e., with singularities at the nuclear sites, 
which tends to zero (or some constant value) at infinity. We are interested in solutions, both scattering (continuum) 
and bound-state, of the Schrodinger equation, 



[V 2 + E - V(r)] i/,(t) = , 
associated with the potential V(r). (As in Ref. [l| we use the Rydberg atomic units.) 



(1) 
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As usual in multiple-scattering theory we shall partition the molecular space as follows: each atomic center i is 
enclosed in a spherical region Tj of radius bi centered on R i7 

7; = {reR 3 : |r-Ri| <h} , (2) 

taken to be mutually disjoint (non-overlapping), so that bi + bj < |Rj — Rj| for all i and j not equal. The surface of 
atomic sphere i will be denoted by ch~i, i.e., 

On = {r el 3 : |r-Ri| = h) . (3) 

We shall consider later a modified partition which includes also an outer sphere r of radius b centered on some 
position Ro, which contains all the atomic spheres, i.e., 

t = {reM. 3 : |r-Ro| < b } , (4) 

such that Tj C ro for all i. In this case the region outside the outer sphere, denoted by tq (the complement of ro), will 
be treated in effect as an atomic center. 



A. The distorted- wave Lippmann-Schwinger equation 

In the usual multiple-scattering methods derivation of the central equations frequently begins with the Lippmann- 
Schwinger equation for scattering states ip( r ) m a multicenter potential V(r), 

yj(r) = Mr) + [ G+(r, s)V(s)^(s) d 3 s , (5) 



where (f>o(v) is a plane wave and G^i^s) is the time- independent free Green function with outgoing- wave boundary 
conditions. (An integral in which the region is not specified should be taken over the whole range of the variable of 
integration; e.g., in the foregoing equation this is therefore over all space.) 

The key to the approach we present here is the splitting of our multicenter molecular potential into two parts: a 
finite non-singular potential Vj(r) which coincides with V(r) in the interstitial region; and a remaining part V^(r) 
which is zero in the interstitial region, but contains atomic potential singularities. Thus, 

V(r) = Vj(r) + V a (t) , (6) 

with 

V IA (r) ifreUti^ 



^f( r ) | V(r) otherwise ^ 

It is not necessary for our immediate purpose to define Via(t) other than that it lead to a non-singular, well behaved, 
and finite Vj(r). One possible choice is to set it to be identically zero. Alternatively, one might choose a form which 
joins continuously or smoothly to V(r) at the atomic-sphere surfaces. Clearly the actual choice will be significant 
when one considers the numerical implementation of the methods we present here and is a point which we shall 
consider in more detail later on. 

Having made this split of the potential our approach to the multiple-scattering problem will be via "distorted- wave" 
theory [l6[, otherwise known as the case of additive potentials, or the two-potential formula [III EH > which we shall 
use here without approximation. In this theory, instead of taking the zeroth-order hamiltonian to be the kinetic 
energy operator, 

H = -V 2 , (8) 
and the scattering potential to be V(r), leading to the usual Lippmann-Schwinger equation ([5]), we take it to be 

Hi = -V 2 + Vj ; (9) 

with the full hamiltonian being 



H = -V 2 +V 



(10) 
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We may then obtain the desired outgoing- wave scattering solutions tJj + of the full hamiltonian by solving the Lippmann- 
Schwinger equation 

ij+(r)= X + (r)+ I G+(v, S )V A (s)ij + (s)d 3 s, (11) 



where G^~ is the Green function relative to the potential Vj, with outgoing-wave boundary conditions, and the 
"distorted wave" x + ( r ) satisfies the Lippmann-Schwinger equation for the "distorting" potential Vi(r), i.e., 

X +(r) = Mr) + J G+(r,s)V7(s) X +(s) d 3 s . (12) 

As is well known (see Ref. [ljj page 132) this may be replaced by the explicit form 

X+(r) = o (r) + / G+(r,s)^(s)^ (s) d 3 s , (13) 



from which, if one knows the explicit form of G^~ , one may evaluate \ + directly. However, our previous work [l| gives 
us precisely this, as we shall see in detail below. 

B. Partial- wave Green-function expansions 

In Ref. [l[ the author has derived an expression for the partial-wave expansions of the exact time-independent 
Green functions, at positive and negative energies, for an arbitrary single-particle potential V(r) (within certain 
limitations). Since the results of this section will be completely general we shall, until further notice, take V(r) not 
to be our multicenter potential, introduced earlier. In fact, the results of Ref. [l|, and some immediate consequences 
which we develop below, will be applied to the more well behaved part Vj(r) of the multicenter potential, providing 
us with explicit forms for Gj (r,s), therefore allowing us to attack the multicenter problem via the integral equations 
(fTTjl and Ipjl. 

The Green functions that we require are solutions of the equation 

[W r 2 + E-V(r)] G(r,s) =5 3 (r-s) , (14) 

together with the appropriate boundary conditions. The potential V(r) is given as a well behaved convergent expansion 
in spherical harmonics about some point which we take as the origin of coordinates, so that 

V(r) =J2Mr)YL(?) , (15) 

L 

where the compound index L represents, as usual, the pair (I, ni) and the summation over L is given by 

oo / 
L 1=0 m=-l 

Again to simplify the algebra we shall make use of real spherical harmonics (where the index m is negative for sine 
and positive or zero for cosine spherical harmonics). As is well known they may be obtained from complex spherical 
harmonics by a unitary transformation which preserves the orthonormality property and the addition theorem. 

The Green functions may be expressed in terms of solutions of the coupled radial Schrodinger equations. These 
may derived from Eq. (fT]) by considering the solution of Schrodinger's equation as a spherical-harmonic expansion, 

^(r) = J2Mr)Y L (r) , (17) 

L 

with which one quickly finds that 

1 ^G' 2 |:))-^ + ^^)-£^'^'W- = ( 18 ) 
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where we have used spherical-harmonic orthonormality and defined the potential matrix elements 

w LL ,{r)=Y,I(L,L\L")v vl {r) (19) 



L" 

(with the Gaunt coefficients for real spherical harmonics given by 

I(L,L',L") = J V L (t)Y L ,(t)Y L „(t)dfi t (20) 

which arc manifestly symmetric under any permutation of the indices L, L' , and L"). It is obvious that the quantities 
wll' are real, for a real potential, and symmetric in the indices. 

The general form of the Green function may be now written as follows |l| : 

G(r,s) = [PL^(r){M- l } L2Li q L , L2 { S )0{ S -r) 

+q LLl (r) {M- l } LiL2 p L , L2 (s)6(r - a)] Y(r) L Y(s) L , , (21) 

where the PLL'{r) constitute a set of real linearly independent vector solutions to the radial Schrodinger equations (| 18|) , 
indexed by V , and regular at r — 0; the qLL'( r ) constitute similarly a set of vector solutions, regular at infinity, whose 
boundary conditions at this limit determine which Green function we have; and Mll 1 is a constant matrix derived 
from the Wronskian of these two sets of functions, as we describe below. Since the potential will, in general, have non- 
zero higher multipole components, the potential matrix wlu will have non-zero off-diagonal elements. Consequently 
the solution matrices will not be diagonal. 

It is interesting to note that the problem of an explicit form for the partial-wave expansions of the Green function 
of a noncentral potential appears to have been rarely treated. One example we have found is Ref. 

[13, 

where the 

authors consider the problem in their Appendix, using the phase-functional formalism with particular boundary 
conditions. They develop an expression without an explicit matrix Ml^i, which seems not readily adapted to the 
further development we shall present later. 

It is useful before proceeding further to introduce some notation. In particular we represent our infinite matrices and 
vectors by sans-serif letters. Therefore we may write equations without explicit reference to indices and summations. 
Thus, in an obvious way, our set of radial matrix solutions may be represented by p(r) and q(r), and the constant 
matrix by M. 

Our general result may be now written more clearly as 

G(r, s) = Y(f) T |p(r) (M" 1 ) 7 q(s) T 9(s -r)+ q(r)M- 1 p(s) T 6»(r - s)l Y(s) , (22) 

where Y(r) is a column vector with elements Y^(r), the superscript T denotes the matrix transpose, and, as was 
shown in Ref. [lj using a well known argument, we have 



p'(r) T q(r) - p(r) T q'M = — M . (23) 



In fact the argument is valid for any pair of matrix solutions, leading to a different constant matrix in each case. 
When the two matrix solutions are the same, further argument shows, at least for p and q individually, that the 
constant matrix is zero, so that 

p(r) T P '(r) = p'W T p(r) , (24) 

and 

q(r) T q'M-q'(r) T q(r), (25) 

provided that, in the latter case, q tends to a diagonal form in the asymptotic limit, which will usually be the case if 
the potential tends to zero quickly enough asymptotically. These two equations would have been automatically true 
in the case of diagonal matrices when the potential is spherically symmetric. 

In the course of the derivation of our main result other relations between these matrices arise. In particular, we 
have 

p(r) (M-!) T q(r) T = ^U^^rf , (26) 
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and 

q'(r)M- 1 p(r) T - p'(r) {M^f q(r) T = 1 1 , (27) 

both of which will be useful later on. 

The choice of particular Green function depends on the boundary conditions, which are imposed on the asymptotic 
behavior of q. For positive energies, E > (with k > such that E — fc 2 ), and Green functions suitable for scattering 
problems, we have two choices for q. Firstly, the choice of outgoing- wave boundary conditions and the retarded Green 
function is obtained with q( + ), defined to be the solutions of the coupled radial Schrodinger equations (TT5|) which 
approach asymptotically, as r — > oo, the corresponding free-electron solutions 

q^ +) (r) = -ih+(kr), (28) 
where the diagonal spherical Hankel-function matrix has elements given by 

h+ L ,(x)=ht(x)5 w . (29) 
We define here similarly spherical Bessel- and Neumann- function matrices j and n, respectively, by 

Oll>{x) = ji{x)5 LL > (30) 

and 

n L L'(x) = ni(x)8 LL , ; (31) 

so that 

h+(x) =j(x) +in(x). (32) 
It is useful to define the diagonal unitary matrices £ and rj by 

[£] w = i l 5 LL > and [ V } LL , = {-l) l 5 LU , (33) 

from which we may see that 

£ 2 = T) and = V 2 = I ■ (34) 

With this notation we may write the famous plane-wave expansion as 

e ik r = 4tt Y(k) T j (Jfer) f Y(r) . (35) 

The choice of incoming-wave boundary conditions and the advanced Green function is obtained with q^ - -*, defined 
to be the complex conjugate of q^ + -*; 

q(-)(r) =q (+) (r)*. (36) 

The constant matrices M from Eq. (|23[1. for the retarded and advanced Green functions, are denoted M+ and Ir- 
respectively. It is easily seen from Eq. (f2"3"]). since p is real, that 

M+ = M_* . (37) 

By the argument which led to Eq. ([23]) we have a similar equation for the Wronskian between q( + ) and q( - \ however, 
the constant matrix may be evaluated by considering the limit as r — » oo, and making use of the well known relation 
for spherical Bessel and Neumann functions [l9j |. 

j[{x)ni{x) -ji(x)n'i(x) = — \ . (38) 



Thus we find that 



q (+)' (r) T q (-) (r) _ q(+)(r) T q (-)' (r) = {3Q) 
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When the potential vanishes we have, of course, the free-electron case and can take p and q to be, respectively, 

p/(r)=j(fcr) (40) 
and (r) from Eq. (|2"8]) . leading to constant matrix 

M +/ = (l/*)l, (41) 

and free retarded Green function 

G+(r,s) = -ikj^ [ji(kr)h+(ks)e(s - r) + hf {kr)j l {ka)6{r - a)] Y L {v)Y L {s) , (42) 

L 

with its famous closed-form representation 

pik\r— s 

G »>' S > = -4^' (43) 

For negative energies, E < (with k > such that E — k 2 ), we have only one acceptable set of functions qj 
regular at infinity, which tends asymptotically to a diagonal exponentially decaying form. In general, the singularity 
of the Green function at a bound-state energy comes about through the singularity of the constant matrix M/ at that 
energy. 

In the free-electron case we can take p and q to be, respectively, 

P/ (r) = i(«r) (44) 

and 

q/ (r) = k+(Kr), (45) 
where, as above, we have defined matrices for the usual modified spherical Bessel and Hankel functions as 

iLL'{x) = ii(x)6 L L> (46) 

and 

k + LL >(x)sk + i(x)S LLf . (47) 

These definitions lead to constant matrix 

M f = -(1/K)r,, (48) 

and Green function 

G (r,s) = [i,(/sr)fc+(«s)0(s - r) + fc+( K r)^( K .s)0(r - a)] Y L (r)Y L (s) , (49) 

L 

with closed-form representation 

To finish this section we note that the standard symmetries of the general Green function (see for example [2(| and 

mi), 

G+(r,s) = G - (s,r)* (51) 

(stemming from the reality of V(r)), and, 

G+(r,s)=G+(s,r), (52) 

may be be verified from the main formula 



8 



1. Further developments 

In the foregoing section we have reproduced the main results from Ref. [l[ together with some useful incidental 
formulas, and set up some notation which will be used in the remainder of this work. In this section we develop 
some immediate consequences of these results. Returning to the positive-energy scattering situation, we derive first 
an important relation between the radial Schrodinger equation solutions p, q^ + \ and q' - -*. 

Since q( + ) and q(~) are the two linearly independent radial-solution matrices regular at infinity, then any other 
solution must be expressible as a linear combination of them. This applies therefore to p. (We note that all three sets 
of solutions are well defined at all radii, determined by inward or outward integration of the coupled radial equations 
from suitable boundary conditions; except for q( + ) and q(~) which may be singular at r — 0.) Thus, given any 
constant vector X (with elements Xl), there must exist constant vectors Z + and ~lr such that, for general r, 

p(r)X = q (+) (r)Z+ + q(-)( r )Z- . (53) 

Evidently we may differentiate both sides with respect to r to get a similar relation between their derivatives. Thus 
we may form the matrix Wronskians of both sides with q^ + -*, q'', and p, successively, and make use of Eqs. (|23jl. 
(|24[). ([25]). and (|39|) . to obtain, respectively, (with some minor rearrangement) three new relations: 

z- = -(|)m + -x, 



Z+ = — M_ J X 



and 



M,Z^ 



.2- = 



(Note that Z + and Z are complex conjugate only if X is real.) We may use the first two of these relations to eliminate 
Z + and Z from our original equation (|53p . and also from the third above, to obtain, noting that X is arbitrary, 



p(r) = 



,(+) 



(r)M. 



>(r)M. 



and 



1+M_ T = M_M+ T 



(54) 



(55) 



From the latter equation we see that the combination M_ M + is symmetric. Since M_ and M + are complex conjugate, 
then M_M + T is hermitian. Therefore it is also real. 

We note, in the positive-energy case, that the relations (j26"|) and (|77|) . between combinations of the radial functions 
and their derivatives, are valid separately for q( + ) and q'~). It will be useful for later reference to prove a further, 
similar relation. So, given p and q, let us define y and z by 



y(r) = rp(r) and z(r) = rq(r) 



Define also u(r), as in Ref. [![, by 



ULL'{r) = w L L'{r) H 5 — o LL , 



ES 



LL> 



noting definition (|19[) and the radial Schrodinger equation (|18p . from which it may be easily seen that 

y (r) = u(r)y(r) and z (r) = u(r)z(r) . 
With these definitions it is not difficult to show, using Eqs. and l|2"7|) . that 



(56) 



(57) 



(58) 



z'{r)M- l )/{r) T - y'(r) (M- 1 ) T z(r) T = I 
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Taking the derivative of this with respect to r, eliminating second derivatives with Eq. (|58[). and using Eq. (|26|) to 
make cancellations, leads to 

z'(r)M-y(r) T =y'(r) (M _1 ) T z'(r) T . 

If we now substitute for y and z from Eq. (|56"|) , and make some cancellations with Eqs. (P6"|) and (ITT)) , then we obtain 
our result: 

q'(r)M-y (r) T = p'(r) (M- 1 ) 7 q'(r) T . (59) 

With our most important result ([54")) from earlier in this section it is possible to eliminate p from the set of relations 
(f!?6")) . (f2T)) . and (|5"9"|) . to obtain a similar set involving only q( + ) and q*- - -*. To carry this out we need first to note that 
Eq. ([55]) may be rearranged as 

M_ T (M^) T = M^M_ , (60) 

showing incidentally that the combination M_ T (M+ 1 ^ T , which will prove important in the next section, is a symmetric 
matrix. One may proceed straightforwardly now to obtain, respectively, 

q(+)( r )q(-)( r ) T = q H ( r ) q (+) ( r ) T , (61) 

q (+)' (r)q (-) (r) T _ q(-)' (r)q (+) (r) T = ( (62) 

which may be compared with the Wronskian (|39p . and 

qW'(r)qH'( r f = q (-)'(r)q(+)'(r) T . (63) 

C. Single scattering and the distorted wave 

Our approach to the multiple-scattering problem is via the distorted- wave Lippmann-Schwinger equation pip . This 
requires, of course, that we have an explicit expression for the distorted wave x + ] which means, in effect, a complete 
solution to the single-scattering problem for the distorting potential Vi. The fact that we have a formula for the 
Green function Gf , from our earlier results, enables precisely that. Since the results of this section can be applied to 
any general, noncentral potential (within the limitations discussed in Ref. [l[), we drop the subscript I, and continue 
to use V to represent an arbitrary potential. 

To find x + then, we use the explicit version (|13[) of the Lippmann-Schwinger equation 

X+(r) - 0„(r) + J G+(r,s)F(s)0 o (s) d 3 s , (64) 

with 

0o (r) = e i;k ' r . (65) 

Now because of the inhomogeneous equation (TTJ} satisfied by G + , and the symmetry we can replace the product 
G + V in the above equation with (V s 2 + E) G + (r, s) — S 3 (r — s). Canceling the inhomogeneous part with the delta- 
function term we therefore obtain 

X + (r) = / 0o (s) ( V s 2 + E) G+ (r , s) d 3 s . (66) 

Since 0o satisfies the free-particle Schrodinger equation, 

(V 2 + E) Mr) =0, (67) 
we can add an extra term to the right-hand side (rhs) to get 

X + (r) = J [0o (s) (V s 2 + E) G+(r, s) - G+(r, s) (V s 2 + E) o (s)] d 3 s . (68) 
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To transform the volume integral into a surface integral, we now wish to use Green's theorem in the form (see 
Ref. i) 



[P(V 2 



E)Q-Q{V 2 +E)P] dr= I (PVQ - QVP) -n do- 
JdT 



(69) 



where dT is the surface enclosing the volume T, and n is the outward-pointing normal to the surface. It should 
be noted here that, since </>o is not a proper normalizable vector in the Hilbert space of physical states, the volume 
integral in Eq. (|68p is not zero. When transformed by Green's theorem to a surface integral, the "surface at infinity" 
provides a non-zero contribution which is the quantity we seek. 

To find this quantity we first consider the volume integral of Eq. as being over a spherical region with large 
radius R, and then allow R to go to infinity. Using Green's theorem then leads to 



sl M/{, 0(s) « [G+(r , s) ]- G+(r , s) ^!)} 



R 2 dQ s 



(70) 



s=R 



where we have made use of the fact that n • V s = d/ds. Now if R is larger than r, then we can take the Green 
function to consist of only the first term of Eq. ([2"2")) . with the choice of outgoing waves q( + ). The plane wave <po may 
also be expressed as the expansion (1351) . Substituting both expressions into the above equation for \ + leads, with 
some rearrangement, to 



X+(r) = 4tt Y(?) T p(r) (MT 1 ) 7, { lim R 2 \q W(i?) T j(fci?) - q^ +) (R) T kj' (kR)] } £ Y(k) , 

I R — >oo L J J 

where we have used the orthonormality of spherical harmonics, 

Jy(s)Y(s) t dfl s = I , 



(71) 



(72) 



to remove the angular integral. To evaluate the limit of the Wronskian expression we make use of the asymptotic 
form of qW and invoke the formula (f38|) . The result for the expression inside the braces is easily found to be (l/k)\. 
Our final result is then 



X+(r) = (47r/fc) Y(f) T p(r) (M; 1 ) £ Y(k) , 
the general form of the scattering state arising from a plane wave in a noncentral potential. 



(73) 



1. Single scattering: complete solution 

Our equation (|73[) contains all the information required for the solution of the single-scattering problem for a general 
potential, but is not usable for this purpose in its present form since it contains the set of radial solutions p regular 
at r = 0. From the scattering point of view we are more interested in the asymptotic region. With this in mind we 
remove p by using the formula f|54[) . This gives us straightforwardly 



X+(r) = 2™ Y(f) T [q<+'(r)A - q (_) (r)l f Y(k) , 

where we define 

A = M_ T (M^) T , 
the symmetric matrix that we encountered earlier. We note therefore that 

At = A* = M+ T (MI 1 ) 7, = A" 1 ; 



(74) 



(75) 



(76) 



i.e., A is unitary. 

Since the asymptotic behavior of x + is more interesting for the single-scattering problem, we consider the form of 
X + at large distances where, by construction, q' + ) and q( - ) approach their free-particle counterparts. Thus, at large 



X+(r) - 2iri Y(r) 



q( +) (r)A 



q^r) £Y(k) 
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If we note that the plane- wave expansion can be written 



- 2m Y(r^ T 

then it is plain that we may write 

X + (r)~0 o (r) + 2^Y(f)' J q^(r) 
We now use the asymptotic form of the spherical Hankel function, 



f q'+'M-qJ-'W £Y(k) 



(77) 



A- I 



£Y(k) 



h+{x) 



as r 



to show that, at large 



We may then write finally 



where we have defined 



_(+) 



(r) 



JkrA 



X+(r) ~0o(r) + 



2?rA 



kr 



Y(r) T [I - A] Y(k) 



'ikr 



A = ^h = e m_ t (m; 1 ) t c, 



(78) 



(79) 



also plainly a unitary matrix. 

Having expressed the asymptotic form of x + m this way we recognize immediately that the usual scattering 
amplitude is given by 



/(k',k) = (27u/fc) Y(k') 7 [l-A]Y(k), 



(80) 



where k' is of the same magnitude as k, but whose direction r is that at which one observes the scattering from the 
incoming wave of direction k. We note here that this result agrees essentially with the corresponding expression in the 
interesting work of Ziegler [22j . It can be fairly easily shown that our matrix A is the same, within some simple matrix 
factors, as his T. His development is a direct approach in terms of asymptotic forms of the solutions of the coupled 
radial Schrodinger equations, rather than, as here, a straightforward application of our explicit general Green-function 
formula. It should be noted, however, that the most important result of our present development is not so much the 
asymptotic form (|78p of x + , which is well known and may be derived from more general principles, but the formula 
(|79| for the matrix A expressing it in terms of directly calculable quantities. 

Of course, it is not always necessary to have an explicit form for A to derive some useful results. It is therefore 
interesting to verify the usual well known properties of / directly from the formula ([8H)) . Thus, to show the reciprocity 
theorem (see [2l[ page 135) we note that 

/(-k,-k') = (2«/fc)Y(-k) T [l-A]Y(-k') 

= (27rVfc)Y(k) T ?7 [I - A] 7?Y(k') 

= (2iri/k)Y(k') T [I - 7 ? A T 77] Y(k) 

= /(k',k); (81) 

where, in the first step we used the parity of spherical harmonics, in the second we took the matrix transpose, and 
we justify the final step by observing that 

V A T V = i^A T ^ v = = A , 

using the symmetry of A and the properties of £ and r\ [Eqs. (|3"3"|) and (|34|) ]. 
Furthermore we have that 



/(k r ,k')7(kr,k) dn kr 



= (47r 2 /fc 2 ) J Y(k') T [I - A f ] Y(k r )Y(k r ) T [I - A] Y(k) dfl kr 

= (4 7 r 2 /fc 2 )Y(k') T [l-A t -A + l]Y(k) 
= (27rz/£0[/(k,kr-/(k',k)] ; 



(82) 
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where, in the first equality we have matrix transposed the (scalar) /* term; for the next step we used spherical- 
harmonic orthonormality, multiplied out the matrix terms and used the unitarity of A; and in the last step split the 
linear sum into two parts. We recognize in this result the generalized optical theorem. As is well known, setting k' to 
k in this result, and observing that the differential cross section is given by 

<7(k r ,k) = |/(k r ,k)| 2 , (83) 

leads directly to an expression for the total elastic scattering cross section, 

a tot = Ja(k r ,k) dQ^ = (An/k) lm/(k,k) , (84) 

the famous optical theorem. 



D. Green-function origin translations and re-expansions 



The mathematical elaboration of multiple-scattering theory relies to a large extent on the translation of expansions 
of various functions, most notably the Green function, from one scattering center to another. Since the Green function 
used is normally the free Green function this poses no problems as it is translationally invariant, and its expansions 
and re-expansions in partial waves are well known. In our case the problem is more complicated since our potential 
need have no translational or rotational symmetry, and this will be reflected in the Green function associated with 
it. It is nevertheless necessary to know how the functional form of the Green function appears seen from different 
centers. It is useful in this context therefore to highlight here a couple of basic mathematical issues that underpin the 
justification for some of the steps in our arguments below. 

The first point that we need to deal with is uniqueness of the Green function. It is, of course, well known that the 
inhomogeneous Schrodinger equation (fT4"]l does not specify the Green function uniquely, since the boundary conditions 
must be imposed. It is easy to show that a solution of Eq. (fT4")l . viewed from a translated origin, will also be a solution 
of the corresponding equation with the original potential viewed from the new origin. Intuitively one suspects that, 
since the boundary conditions are imposed at r — > oo, and are the same in both cases, the resulting solutions will be 
just different points of view of one unique Green function; however, this needs to be shown formally. 

In fact it is the Lippmann-Schwinger equation for the Green function, 



or equivalently, 



G+(r,s)=G+(r,s) + J G+(r,t)V(t)G+(t,s) dh, 
G+(r,s)=G+(r,s) + J G+(r,t)F(t)G+(t,s) dh, 



(85) 



(86) 



which guarantees uniqueness (see Ref. [18| page 133 or for a detailed proof Ref. [20] page 490 ff). 

It is not difficult to check that the unique solution G + of Eqs. ([83)1 and (1831) . when considered relative to a translated 
origin of coordinates, satisfies the corresponding equation in the new coordinate system. (Of course, the potential 
V will also have a different functional form in the new coordinates.) This uniqueness guarantees that when we find 
a form for G + in some coordinate system, it will be equal to a translated version of the form arising with another 
origin. 

It is useful to check directly that our general form (f2"2"]l does indeed satisfy Eqs. ([85| and (|86[) . Thus let us assume 
G + is of the form of Eq. ([22]) and consider Eq. ([85)) for example. Then we should have 

G+(r,s) = G+(r,s)+ J G+ (r,t)V(t)G+ (t,s) dh 

= J G+(r, t) (V t 2 + E) G+(t, s) dh 

= G+(r,s)+ J [G+(r,t)(V t 2 +£)G+(t,s)-G+(t,s)(V t 2 +£)G+(r,t^ dh 



= G+(r,s)+ 1,1m J |G+(r,t)| [G+(t, s)] - G+(t, s)| [G+(r,t)] 



R 2 (Kit; 



t=R 
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where, in the first step we have replaced the VG + in a similar way to that of Eq. (|66f [which is justified because our 
form for G + was developed directly from Eq. (|14p] , in the second we use a similar property of Gj , and in the third 
we invoke Green's theorem for the region within a sphere of radius R (taken to be greater than r and s) which we 
allow to go to infinity. Clearly, for the equation to be true we require that the limit of the integral be zero. Using 
the appropriate expansions for the Green functions, noting that t = R > r, s, we may rewrite this (making use also of 
spherical-harmonic orthonormality) as the expression 

M+- 1 p(s) T Y(s) . (87) 

Now, since qW was chosen to be asymptotic to the free-particle case qy + ' (assuming as usual a potential which tends 
quickly enough to zero), the Wronskian will tend to zero as desired. So, with the foregoing arguments in mind, we 
may have confidence that the form of G + that we derive in some local coordinate system is identical to that found in 
some translated coordinates. 



1. Re-expansion of free Green functions 

Let us now derive the re-expansion formulas for the free Green function Gq . These are, of course, well known (and 
may be found in Ref. [2| whose approach we follow closely). However, it is useful to repeat briefly the derivation, 
since it allows us to introduce some new notation and express the formulas in a different form than usual. Also it is 
a model for our later derivations. 

First, we recall the atomic centers Ri (for i = 1, . . . , N) which we introduced earlier. For some general position 
vector r let us denote the relative position vector with respect some given atomic center Ri by 

rj = r — Rj . (88) 

For the purpose of developing the re-expansion formulas we may consider, say, two distinct centers, Ri and Rj, as 
being two general points in space. The resulting formulas will then be immediately useful for our multiple-scattering 
problem. 

We consider first the re-expansion of the free radial Schrodinger-equation solutions, p/(r) = }(kr), and begin by 
considering a plane wave expressed in terms of coordinates relative to Ri and Rj thus: 

exp(ik • r) = exp [ik • (r, + Ri)] — exp(ik • r^) exp(ik • Rj) 

and 

exp(ik • r) — exp(ik • Yj) exp(ik • Rj) . 

Equating the two expressions and using the plane-wave expansion (|35|) leads, with some rearrangement, to 

Y(? t ) T i(k n ) £ Y(k) = YfofXfcr,) £ Y(k) exp(»k ■ Ry) , 

where we have defined the atomic center-to-center vector by 

Ry = Rj — Ri . (89) 

If we multiply both sides of our equation on the right by Y(k) T and integrate over all k, then we get 

Y(r l ) T j(fcr 4 ) = y(?j) T ](k rj ) D(k; Ry) , (90) 

where we have used spherical-harmonic orthonormality and defined the matrix associated with displacement, for 
general R, by 

D(fc;R) = ffi(A;R)ft (91) 

and 

D(jfc; R) = f Y(k) exp(ik • R)Y(k) T dn k . (92) 



Y(rf P/ ( 



M 



lira 

R—*OG 



U;\Rf q ^'(R) 



q ( f +) (i?) T q (+) 



R 2 
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If we use the plane-wave expansion for the exponential in the last equation, evaluate the angular integrals, and 
rearrange, we find that the matrix elements of D(/c; R) are given by 

D LL , (k; R)=4ttJ2 (kR)Y L „ (R)I(L, L' , L") , (93) 

L" 

from which one may note that D(fc; R) is symmetric. The matrix elements of D(fc; R) are now easily seen to be 

D LL ,(k;R) =4irY / i l "~ l ' +l Ji"(kR)Y L , / (R)I(L,L',L") . (94) 

L" 

If we now note that the Gaunt coefficients I(L, L' , L") are non-zero only when I + 1' + I" is an even integer, then we 
see that D(fc;R) is real. Writing out fully the summations in Eq. (|9"0|) , and using the explicit matrix elements (I94p . 
allows us to recover the usual re-expansion formula ■ 

Now since the points indexed i and j are arbitrary we may interchange them to obtain another equation like Eq. 
(|90|). which we use to eliminate the term Y(r J ) T j(fcr : ,). On rearrangement we then find 

YfofUkn) [I - D(fc; -Rij)D(k; R y -)] - . 

Multiplying on the left by Y(?i), integrating over all angles, and invoking spherical-harmonic orthonormality, then 
leads to 

j(fcr-i) [I - D(k; -R,j)D(fc; %)] = . 

Since the diagonal matrix \(kri) consists of the linear independent (regular) solutions of the radial Schrodinger equa- 
tion, then the expression in square brackets must be the zero matrix. In other words, since R^ is arbitrary, then we 
have that, in general, 

□ (fcjR)- 1 = D(fc;-R) . (95) 

Furthermore, we have 

D^R)- 1 = £,D(k;-R)^ 
= ZD(k;Ti)*? 
= £D(fe;R)t^ 
= D(fc;R)t , 

where, in the first step we take the complex conjugate of Eq. (|92[) . and in the second we invoked the symmetry of 
D(k; R). Thus D(fc; R) is unitary. Since it is real, then it is orthogonal. It may be also easily seen that 



D^R)- 1 = D(fc;-R) = D(/c;R)t . (96) 

(+ 
/ 



To develop the corresponding formula for the free (outgoing- wave) radial solutions regular at infinity, qi + ^(r 



— ih + (fcr), we need to use the translation properties of Gq . However, since these solutions are singular at the origin, 
the re-expansion must be approached differently above and below the separation distance, which regions we shall 
denote as the "far" and "near" regions, respectively. 

With reference to the diagram in Fig. [TJ let us consider a general point with position vector s, on the surface of 
a sphere centered on Rj, with a radius somewhat greater than Rij. Let r represent a general point on a sphere of 
positive radius centered on R^, contained entirely within the first sphere. Such a sphere can always be found, since 
Sj > . It may be then seen that in this configuration we shall always have Si > ri and sj > Tj . Therefore we may 
write 

G+(r,s) = - f^' 5 '' = -ifcY(f l ) T j(fcn)h + (fc Si )Y(s l ) . 

47T \Ti - 8i\ 

Translation invariance allows to write also, since — s.^ = — Sj, that 
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Equating the two right-hand expressions gives us 

Y^) T j(*ri)h + (*«i)Y(Si) = Y(r j ) T j(kr j )h+(ks j )Y(s j ) , 

which may be rewritten, using Eq. (|90p . as 

Y(?i) T j(*Ti) [h+(fc Si )Y(s 2 ) - D(fc;-R y )h+(fc Sj )Y(%)] = . 

Since the sphere associated with r is entirely within that associated with s, then the last equation is true for all r.j 
and at least a small range of r». So we may invoke our previous argument to conclude that the expression in square 
brackets (which has no dependence on r^) is zero. Thus, transposing and using the properties of D(/c;R), we may 
write 

Y(s i ) T h + (fcs l ) = Y(s j ) T h + (ks j )D(k;R ij ) for s 3 > R tj , (97) 

the far-region re-expansion formula. This may be compared with Eq. (|90p (indeed the real part corresponds identi- 
cally) . 

To obtain the re-expansion formula for the near region we first write out Eq. (|97[) explicitly, so that we have 

h+(ks t )Y L {s t ) = 4tt ^(kR^h+iks^YL^R^Y^is^IiL, L' , L'') ; (98) 

L'L" 

where s, = Sj + R^- with the restriction Sj > Rij. In the near-region case, where Sj < Rij, we may use the same 
formula with the variables interchanged, since the vectors are arbitrary apart from the order relation. Therefore, we 
have 

h+(ksi)Y L (si) = 4tt il "~ l+V ' 3v'{ks j )hi l {kR ij )Y L , l {s j )Y L ,{%)I{L,L l ,L n ) . (99) 

L'L" 

This may be written in our notation as 

Y(s,) T h+(fc Si ) = Y(s J ) T j(A; Sj )F(fc;R y ) for 8j < R. l3 , (100) 
where we have defined F(fc; Ry) by 

FLL'(k; Rjj) = ATTj2^ l ' +l "K,(kR l j)Y L „{n i3 )I(L,L\L") . (101) 
V 

2. Gaunt algebra 

Having encountered the Gaunt coefficients yet again in our formulas it seems worthwhile to investigate them in 
more detail at this point. We shall look briefly at some of their properties, but consider them from a slightly unusual 
point of view. 

Normally the Gaunt coefficients arise in the re-expansion of the product of two spherical harmonics, where 

Y L {v)Y L ,{v) =Y J I{L,L',L")Y L „{v) . (102) 

L" 

Use of spherical-harmonic orthonormality leads immediately to the definition (|20[) that we presented earlier. In the 
case of complex spherical harmonics the corresponding versions of the coefficients have a well known expression in 
terms of Clebsch-Gordan coefficients. 

The particular form of the definitions of matrices such as w of Eq. (|T9|) . D(fc; R) of Eq. (|92|) . and F(k; R) of Eq. 
(|100p . for example, suggests that it may be useful to consider a set of "Gaunt matrices" f~( L \ indexed by L, with 
matrix elements defined by 

r( L '\ L , =I(L,L',L") , (103) 
easily seen to be real and symmetric. We may now reinterpret Eq. I|102p as an eigenvalue equation, 



[T L (r)]Y(r) = T^YC?) , 



(104) 
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in which, for any r, the vector Y(r) is an eigenvector of the matrix f~( L \ with eigenvalue Yz,(r). 

Now the properties of the Gaunt coefficients allow us to say something about the structure of the individual matrices. 
In particular, we note that the triangle condition on I, I', and I" implies that the Lth row (column) of V^ L ) can only 
have a finite number of nonzero entries, from index V — \l" — l\ up to index V = I" + I. Furthermore we note that the 
matrix is banded with respect to I. 

Of special note is the fact that the matrix elements of the first row (column) are given by 

r< L "\o, 0) L> = I ((0, 0), L', L") = (47r)- 1 /2 (5i , i „ , (105) 

from which we may deduce that the Gaunt matrices are linearly independent. 

For a given vector a with (in general complex) components a^, let us define a linear combination of the Gaunt 
matrices by 

A[a] =^2 a L r(L) ■ ( 106 ) 

L 

There is no need to consider convergence in this definition since each matrix element of A [a] only contains contributions 
from a finite number of the ol, by virtue of what was said above. Linear independence gives us that 

A [a] = a = . (107) 

If one is given a A [a], one may recover a by reading off the first row (column) and using Eq. (|105p : which shows that 
the correspondence between a and A[a] is one-to-one. It is also not difficult to show that A[a] is a linear function of 
its parameter a, and that, if the products exist and are finite, 

A[a]b = A[b]a. (108) 

Among the cases that we have encountered we may write, for example, 

D(fc; R) = A[4<j(fci?)Y(R)] . (109) 

To see some more interesting properties let us consider now the integral of the product of four spherical harmonics, 

h(L, L', L", L'") ee J Y L (t)Y L , (t)y £ „ (t)Y £ ,„ (t) dfi* , (110) 

which are manifestly symmetric under all permutations of the indices. We may use Eq. f|102[) to substitute for pairs 
of spherical harmonics in this definition, then use orthonormality to obtain three possible equivalent expressions, so 
that 

h{L,L',L",L"') = 2/(L,i" / J i 1 )/(L' J L" J L 1 ) 

Li 

= Y i I{L,L",L 1 )I{L"',L',L l ) 

Li 

= ^2 L' , Li)I(L" , L'" , Li) . (Ill) 

Li 

Making a particular choice for the assignment of indices leads to the matrix equations 

V {L) V (L>) = V {L') v iL) = I( L , V, L,)^ . (112) 
Li 

In summary then, the form a set of real, symmetric, linearly independent, mutually commuting matrices, the 
product of any pair of which may be expressed as a linear combination of some of them (in fact in the same way as 
the spherical harmonics from which they derive) . The fact that they commute with each other is, of course, consistent 
with their having common eigenvectors. 

It is interesting now to consider the r( L ) as generators of some larger set of matrices. In particular, let us denote 
by Aq the set of all linear combinations A[a], i.e., the set spanned by the r( L ). By virtue of the relations (| 1 1 2[) we see 
that this set is closed under matrix multiplication; although, when we consider the product of two elements of Aq, 
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we must be careful here about convergence. Indeed, we see that the general matrix element of a product of linear 
combinations may be written 

which, to be finite, requires that the admissible vectors for Aq represent convergent spherical harmonic expansions of, 
say, square-integrable functions of f . Clearly Aq is a commutative algebra of matrices which we shall call the "Gaunt 
algebra". (We conjecture also that is possible that the I = 1 Gaunt matrices generate all the rest.) 
It is not difficult to show, using the equations (jllljl . that we have 

A[A[a]b] = A[a]A[b] . (113) 

This may be used to show a further property of D(fc; R). To see this rewrite the re-expansion formula (|90D as 

eXArOYft) = B(fc;R iJ 0£j(fcr J -)Y(r j ) . 

Using Eq. (|109j) we see straightforwardly that 

D(fc;r,) - A[4<j(fcr 4 )Y(f,)] 

= A[D(fc; Ry)47r£j (krj )¥(?,•)] 
= D(fc;R 4 ,)A[4<j(fcr,)Yfc)] 
= 6(jfc;Ry)D(fc;rj) . 

Now Yi = JUij + Yj with no restrictions, and the vectors are effectively arbitrary. So we see that, for general vectors 
R and S, we have 

D(/fc;R + S) = D(/c;R)D(fc;S) , (114) 

and it is easy to see that the corresponding relation holds for D(fc; R). We note here that this relation can be proved 
directly from the definition (|92p if one makes use of the angular delta-function property of spherical harmonics, 

J2Yl(t)Y l (s) = Y(r) T Y(s) = 6(r -s) . (115) 

L 

This property has been noted before by Natoli et al. in Ref. [23| [see their Appendix A noting that our D^ik] Ry) 
corresponds to J^ L , in the notation of their Eq. (A23)]. 

3. Re-expansion of general regular solutions 

We may now derive the re-expansions of the general coupled radial Schrodinger-equation solutions p and q. For the 
general regular solutions p we consider the explicit form (|13[) of the Lippmann-Schwinger equation and the resulting 
scattering state x + ( r ) 01 Eq. (jT3"l) . noting that the plane- wave term <po(r) occurs linearly on the rhs of the former. If 
we now consider an initial plane wave fao (r) given by 

&o(r) = e lk - r * = e 4k - r e- 4k - R * = e - ik ' R * </> (r) , (116) 

then we may repeat the derivation of Eq. ([73")) with respect to a new origin at R,;. This leads to a scattering state 

Xt(ri) = (47r/fc) Y(f l ) T p,(r l ) (M7 + 1 ) T fY(k) , (117) 

where the subscript i on p^, Mi + (and implicitly q^) shows that the coupled radial Schrodinger-equation solutions 
with respect to the center R^ will be, in general, different to those of other centers. If we now invoke the linearity 
with respect to (pqM of the rhs of Eq. (|13[) . and the uniqueness of solutions of the Lippmann-Schwinger equation (see 
for example Ref. [20(), then we may conclude that 

X + (r) - e^xtin) . 



(A[a]A[b]) LL , = / Y L (t) 



L'" 



b L ,„Y L „,(t) 
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(Since the initial plane waves differ only by a constant phase factor, it should be intuitively obvious that the scattered 
waves must be essentially the same within a phase factor.) Therefore we have 

X + (r) = (4*/*) Y(r,) T p^) (M- 1 ) T CY(k) e lkR * . 

A similar result may be derived with respect to some other center Rj, and we can equate the two expressions to 
obtain 

Y(r l ) T Pi (n) (Mr + i) T £Y(k) e *'*< = YftfPifo) (M^f £Y(k) e ik ^ . 

From this point we may proceed as in the free Green-function case and find that 

Y(r,) T P ,(r l ) (Mr + i) T = Y(r i ) T p i (r i ) (Mj^f D(*;Ry) , (118) 

with no restrictions on the vectors involved. This is then the re-expansion formula for the general regular solutions. 
It is easy to see that it reduces immediately to the result (|50")l when then scattering potential vanishes. 

We make two observations about our formula (|1 18[) . First, it may be noted that the radial solution matrices p are 
required only to be linearly independent, real, and regular at the origin. We do not specify any further boundary 
conditions or particular form at r = Q. At first sight this may be thought to pose a problem when we relate 
solutions centered on different origins whose boundary conditions may be arbitrarily chosen. In fact this problem is 
avoided because of the presence of the inverse of the constant Wronskian matrix M for each center, which serves to 
"standardize" the regular solutions by reference, through the q ( - + - ) , to the common boundary conditions at infinity. 

Second, we note for the free radial functions that the re-expansion formula has a certain symmetry because both 
the functions themselves and the displacement operator D(fc;Rjj) involve the spherical Bessel functions. In fact this 
symmetry underlies the product formula (| 1 14|) . As a consequence the expressions on both sides of the formula may 
be easily seen to satisfy the free Schrodinger equation when considered as functions of r j , rj or Ry . If one wishes to 
investigate the analogous properties of our latest expression, then one must take into account the "hidden" dependence 
of Pi on Ri. 



4- Re-expansion of general irregular solutions in the far region 



With the result of the previous section we may proceed quickly to obtain the re-expansion formula for general 
irregular solutions in the far region. To see this we consider again Fig. [TJ where the position vector s is such that 
Sj > Rij , and we have Sj > and Sj > rj . Expressing the general Green function in terms of vectors relative to the 
centers Ri and Rj in turn, using Eq. (|22[) . leads to 



G+(r,s) = Y(f l ) T p l (n)(Mr + 1 ) T q l (+) ( S T Y(s l ) 

= Yi,v' P ,«r,i (MT + i) T q (+)( Sj ) T Y(s J ) . 
Making use of the previous section's results leads immediately to 

Y(f J ) T P,(r l ) (M^f [q«( Si f Y&) - D(fc; -RyJq^tofYfo 



(119) 



= 



By similar reasoning to that which led to Eq. (f97|) , we see that the expression inside the square brackets in the last 
equation must be zero, and we may then write (on transposing) 



Y(s i ) T ql +) (s i )=Y(s j )' r q)- rj (s j )D(k;R ij ) for Sj > R t 



\Tj+) 



(120) 



the general far-region re-expansion formula. It is also easily seen to reduce to the free result when the potential 
vanishes. Taking the complex conjugate of (|120[) . noting the reality of the spherical harmonics and the D(fc;Rjj), 
leads to the corresponding formula for the incoming-wave solutions: 



Y(s i ) T q(- ) ( S ,)=Y(%) T q(-)( Si )D(fc;R ij -) for Sj > R l3 



(121) 
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5. Translation of the single-scattering amplitude 



Our re-expansion formulas for the general irregular solutions in the far region allow us to take the development of 
the single-scattering problem further. In particular we may now elucidate the behavior of the scattering amplitude as 
seen with respect to different origins. Intuitively, since the scattering is calculated from the asymptotic behavior of 
the solutions, one expects that a finite displacement should have no real effect. This is not quite the case as we shall 
see. 

First, however, we must use the formula (|54p to rewrite the general Green function (|22j) in terms of only the irregular 
solutions. With some algebra we find, noting (|75|) and (p0|) . that 

G(r,s) = 

yY(r) T [q( + »(r)Aq( + '( S ) T - q(~) (r)q<+> (s) T 9{s - r) - q(+\r)q^(s) T 6(r - s)] Y(s) . 

(122) 

We now consider this form expressed relative to the centers and Rj in turn, as before, to give us two equivalent 
forms for G(r, s). Let us then take the general vectors r and s to be such that rj > i?y and Sj > Rij, and use the 

re-expansion formulas (|120[) and (|121[) to replace the qj + ^ and q terms in the first of these equations. We should 
also choose r and s to be on spheres centered on Rj such that, say, rj < Sj by a margin large enough so that r, ; < Si 
as well. This ensures that the theta functions correspond. Equating the right-hand sides, canceling the theta-function 
terms, and rearranging, we find that 

Y&fq^fo) [D(A;R«)A i D(A; -Ry) - A,] qj+'fefY^) = . 

For large enough rj and Sj their order relations above are true for all values of the angular variables, so we can 
conclude that the expression in square brackets is zero. Using our earlier definitions we then find that 

D(k; Ry )A l D(fc; -Ry) = Aj . (123) 

We have seen in Eq. (|109[) that D(fc; Ry) is a linear combination of Gaunt matrices and as such shares the common 
eigenvectors Y(r) for arbitrary?. It is not difficult to see that the associated eigenvalue is exp(i/cr • Ry), and so 

Y(k') T A,Y(k) = Y(k') T D(fc;Ry)A l D(fc;-R l3 )Y(k) 

= exp[i(k' - k) • Rij] Y(k') T A,Y(k) (124) 

If we now consider the term from Eq. (|80|) for the scattering amplitude, involving the identity matrix, we see that it 
is precisely of the form of Eq. (|115[) . i.e., a delta function of the angular variables. Therefore it is unchanged if we 
multiply by exp[i(k' — k) • Ry], since this is a function of the difference of the same variables and is unity when they 
coincide. So finally we may write 

f 3 (k', k) = exp[i(k' - k) • Ry]/;(k', k) . (125) 

We see from this that the cross section is unchanged, but there is a phase shift when the scattering amplitude is 
calculated with respect to a different origin. The first is entirely in accord with our expectations, while the second is 
what gives rise to interference effects if we had several disjoint scattering potentials. 



6. Re-expansion of general irregular solutions in the near region 



In the case of the near-region re-expansion formula, we do not yet have a generally applicable method leading to 
a closed formula similar to those of the immediately preceding sections. The argument leading to the re-expansion 
formula (|100[) for the free Green function is not available to us in the case of a general potential, since it makes use 
of the fact that the free radial solutions are independent of which center they are referred to. 

Let us therefore look for an alternative approach by considering the situation in Fig. [51 returning to the free Green 
function in the first instance. Thus we choose an arbitrary general point s on some sphere centered on Rj with, now, 
sj < Rij; and an arbitrary r on a sphere centered on R^, entirely outside the first sphere. Since here we shall always 
have Si > ri and Sj < r.j , an argument similar to that preceding Eq. (|97[) leads instead to 

Y(?. i ) T j(fcr l )h+(fc Sl )Y(s. i ) = Y(?,) T h+(fcr 3 )j(fc S ,)Y(s,) . (126) 
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One sees here already, since r and s are independent and arbitrary, the likely proportionality between h + (ksi)Y(si) 
and ](kSj)Y(Sj). To see this explicitly we first multiply on the left by Y(?j), integrate over the angular variables fj, 
and use orthonormality of spherical harmonics to get 



j(fcr i )h+(fcs i )Y(s i ) = 



Y(r i )Y(r i ) T h+(fcr i ) dtt ri 



}(ks 3 )Y(s 3 ) 



(127) 



To eliminate the }(kri) term, and thus the dependence on on the left-hand side (lhs), we use the famous Wronskian 
formula (|38|) for spherical Bessel functions in the form 



h+(kr)-^ [j(fcr)] - ^ [h+(fcr)] j(fcr) = -^1 



(128) 



Thus, if we multiply our main equation on the left by the derivative with respect to (wrt) of h + (fcri), and the same 
derivative of our main equation by h + (fcri), and take the difference, we obtain, after some rearrangement, 



h + (fcs,)Y(si) = ikK^ T j(ksj)Y(sj) for Sj < i?y ; 



where we have defined the important matrix 



K 



J* T = 2 



h+(fcr,)Y(f J )^- [Y(r 3 ) T h+(fcr J )] - A [h+(fcr,)] Y(f 4 )Y(r,) T h+(fcr 3 ) }> dO ri 



(129) 



(130) 



If we introduce some notation from Ref. [23j for vectors of free radial solutions, denoted J(fcr), N(fcr), and H ± (fcr), 
with elements given by 



Ji(fcr) = ji{kr)Y L {v) 



(131) 



and 



then we may write 



N L (kr) = ni(kr)Y L (r) , 



H±(kr)=hf{kr)Y L {r) 



(132) 
(133) 



dn (ri ) 



[H+(fcr,)] H+(fcr,) T - H+(fcr,) A [H+(fcr,) T ] J dft ri 
{V rs [H+(fcr,-)] H + (fcr,) T - H+(fer,-)V r< [H+(fc ri ) T ] } • n da 



_0_ 

dr 



(134) 



where the volume Ti(r{) is a sphere of radius centered on R^. 

Now Eq. (|129p is clearly another form of the near-region re-expansion formula (jlOOp for the free irregular solutions, 
and comparison with this latter shows, using also Eq. (|106[) . that 



= -(i/k)F(k;Hij) = A[4<H+(fcR y )] ^ 



(135) 



It should also be clear that Kj l is constant, even though the definition (|130|) and alternative forms (|134l) appear to 
depend on r$. We may show this directly by using Green's theorem to transform the second form in Eq. (|134[) into 
a volume integral. First, however, we note that the gradient operator in the integrand of this latter form may be 
taken to be wrt r, since is just a translation of this by a constant vector. One then sees that this integrand is 
anti-symmetric under the simultaneous exchange of the indices i and j and matrix transposition. Second, from Fig. [2] 
we see that the constraints on our geometric configuration imply that r» < Rij = Rji and we can therefore use our 
preceding results (|129p to write 



h+(fcr J -)Y(f J ) = ik\C] T j(Ar<)Y(ri) 



(136) 



We may now use this and Eq. (|129p to eliminate the irregular solutions altogether from Eq. (|126p and find, on 
rearranging, that 



Y(r,) T j(fcr l ) 



ikKi - ikK{ 1 T 



j(fc S ,)Y(s J ) = 
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from which we deduce that 



Now the integrand of Eq. (|134|) is singular at r = Rj and r = Rj . Using Green's theorem 



(137) 



we find that 



K = i { [( y2 + E ) H+ ( kr j)] H + (fcr,) T - H+ikrj) [(V 2 + £)H+(fc ri ) T ] } d 3 r . 

Since the vectors H + (fci"j) and H + (fcrj) are solutions of the free Schrodinger equation, the integrand is identically zero 
away from the singularities. Thus the only finite contribution to the integral comes from the singularities themselves. 
It is then clear that the surface integral expression for K^ 1 may be taken over any surface which encloses only the 
singularity at r = Rj. In particular, the integral is independent of r*j. The possibility then arises that the closed form 
of Eq. (|135|) may be obtained by evaluating the integral as — > 0. This is not so interesting in the free-electron case, 
but might be applicable to our general irregular solutions which we consider below. 

It is not difficult to show that the contribution to such an integral from the singularity at r = Rj would be equal 
and opposite to that at r = Rj . This may be seen by evaluating the second form in Eq. (] 1 34[) over the surface drj (rj ) , 
using the anti-symmetry of the integrand to show that the result is — T , and noting Eq. (| 1 3T[) . Alternatively we 
might consider an integral over a surface which encloses both singularities and show that the result is zero. One simple 
possibility is a spherical surface of radius R > R^ centered on r = Rj, so that we can use the far-region re-expansion 
formula (j9"Tj) to replace the term H + (kYj). Thus we would have, say, 



[ {V [H + (fcr 3 -)] H + (fcr,) T - H + (fcrj)V [H+(/crj) T ] } ■ n da 
= D(A;Ry) / {V[H+(fc ri )] H+(fcr i ) T -H+(fer i )V[H+(fer i ) T ]}-nda 

Jdn(R) 

= 0, 

where the last step is justified by noting that the integrand is anti-symmetric and, when the integration over angles 
is done, the result is diagonal. 

With the foregoing development in mind we may proceed to consider in a similar way the near-region re-expansion 
of the general irregular solutions. Noting again the situation in Fig. [5] and the order relations, the expressions for the 
general Green function relative to the two centers, Rj and Rj, lead to 



(138) 



Y(f i )Y(f i ) T qf ) (r i ) dQ T 



As before we have then 

P,(r J )(Mj + 1 ) T q| +) ( Sl ) T Y(s l ) = 

Using now Eq. I|23p we obtain 

q| +) ( Sl ) T Y(s,) = K- 7 '* T (MT + !) Pj ( Sj ) T Y&) for Sj < R ij ; 

with the definition 



(139) 



qW(r 4 ) T Y(f l )Y(? J ) r A 



q5 +) (r,-) 



_d_ 



(140) 



One may easily check the correspondence with our earlier formula in the limit of zero potential. Again we may consider 
the re-expansion from center i to center j and obtain 



q$ +) (r i ) r Y(f i ) = K« T (M.^ 1 ) p. i (r 4 ) T Y(f l ) for n < ft. 
using the two re-expansion formulas to prove the symmetry 



(141) 
(142) 
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Let us define vectors of general solutions by 

P(r) = p(r) T Y(f) , (143) 

and 

Q (±) (r) EEq( ± )(r) T Y(r) . (144) 

We may write 



K ]l = -r, 2 



f \v\Q^\v j )]Q^\v i f-Q^\r j )w\Q\ + \r i f]}.nda, (145) 



and show, as above, the independence wrt r^. Precisely the same comments as before, about surface and volume 
integrals, apply in this case, except that, in using Green's theorem it should be noted that the terms like (V 2 + 
E)Qj(Tj) are not zero, but Vj(Tj)Qj + \rj), and they will cancel each other out in the integrand of the volume 
integral (since Vi(ri) = Vj(rj) = V(r) by definition), so that this is again zero away from the singularities. 

In the absence of a closed formula for the K Jl we must evaluate them numerically via Eq. (| 145|) . The freedom 
in choice of surface may be invoked to use the atomic-sphere surfaces of radius bi that we introduced earlier. The 
possibility that a closed form might be obtained by taking the limit as — > remains to be investigated. 

E. Multiple-scattered-wave method 

Having developed the necessary new mathematical machinery we may now return to consideration of the multiple- 
scattering problem. Here we follow the general approach of Ref. but with several important differences. It will 
be clear that the foregoing development of re-expansion formulas for general solutions of Schrodinger's equation is 
intended for the distorting potential Vj(r). Thus it is assumed that we have available, for each atomic center with 
position vector Hi (for i = 1, • • ■ , AT), the exact solution matrices pi and , regular at Vi = and oo, respectively, 
in this potential. Finding these solutions is thus, of course, a large part of the computational burden in our approach 
to the multiple-scattering problem. For each center we evaluate the constant Wronskian matrix M; + from these 
solutions, via Eq. (|23[) . at some suitable radius ri = bi, say. 

We are interested in finding the scattering solution (r) of the Schrodinger equation for our multicenter potential 
V(r), which develops from an incoming plane wave </>o( r )- It is assumed that we have solved the scattering problem 

for our potential Vj(r) in terms of the pi, q[ , and Mj + just mentioned, and thus have expressions for the Green 
function G/(r,s), of the form of Eq. (|22[) , with respect to any atomic center i. In particular, we have an explicit 
expression for the distorted wave x + ( r ) given by Eq. ([73)) 

Following our earlier discussion it is then clear that to find ip + ( r ) we must solve the Lippmann-Schwinger equation 
(fTTjl with the Green function G/(r,s) and the singular part Va(v) of our multicenter potential, together with the 
distorted wave as inhomogeneous part, i.e., 

^+(r) =x + (r) + J G+(r,s)Vk(s)V + (s) d 3 s . (146) 

As usual the first step is to split the volume integral in this equation according to the molecular partition. Therefore 
we write 



V^+(r) = X +(r) + 
where we have denoted the interstitial volume by I, so that 



V ( G+(r,s)y A (s)^+(s) d 3 s+ ( G+(r, S )V A (s)i> + (s) d 3 s , 
i=i J n Jx 



N 



(147) 



Now, by construction, the singular part of the potential is zero in the interstitial region, so the last term in our 
equation vanishes, and we therefore have 



N 

iP + (r) = X + (r) + J2 / G}(r lS )V A (s)iP+(s)d 3 S . 
i=i Jt - 
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Now ip + (r) is a solution of the Schrodinger equation |T]) for the full potential, and because of Eq. (J6]) we have 

VA(r)V+(r) = [V 2 + E - Vr(r)] ^+(r) . (148) 
Using this we can eliminate Va(t) and obtain 

N 

V+(r) = X+W+E / G+(r,s)[(V s 2 +i : ;)^+(s)-U / (s)^ + (s)] d 3 s 

- X + (r)+]T / {[G+(r,s)(V 2 +i^ + (s) 

-Tp + (s) (V 2 + B) G+(r, s)] + 5 3 {r - s)^+(s)} d 3 s , (149) 

where in the last step we have used the basic property (fT4|) of the Green function, noting the symmetry (|52|) . Because 
of the delta-function term we see that there are two distinct cases to be considered. Invoking first Green's theorem, 
we find that 

V+(r) = X + (r)+^ / [G+(r,s)V s V> + (s)-^ + (s)V s G+(r,s)] .^do % ifr^(Jr,, 

%=\ Jdr - i=i 



(150) 



and 



N N 

= X + (r)+]T/ [G+(r,s)V s ^ + (s)-V^ + (s)V s G+(r,s)] ■n l dcj l ifre|J r - 



(151) 



We may now rewrite these equations in terms of radial Schrodinger-equation solutions at each atomic center. We 
have already discussed the Green function at considerable length, and the associated radial functions may be taken 
as known. We require also expressions for the wavefunction ip + ( r ) inside each atomic sphere r^. Here we follow the 
development of Natoli, Benfatto, and Doniach [5j who show that we may write, for r 6 Tj, 



(152) 



LL' 



where the C l L , are constants, and the R^t'^i) constitute a set, indexed by L', of linearly independent solutions to 
the radial Schrodinger equations at center i, 



1 / d 



d 



T7 \ dri V 1 dr 



1(1 + 1) 



E 



which are regular at the origin. The potential matrix is given by 

vi L ,(n) = J2i(L,L>,L")vi„(n) 



(153) 



(154) 



where V^{ri) comes from the partial-wave expansion of the full molecular potential V(r) at center i, i.e., 

l/(r)^^V2(r 4 )r L (? 2 ). (155) 



In our matrix notation we may write the atomic solution as 

</>+(r) =Y(?,) T R l (r 4 )C 



(156) 



If we now consider the case where our general vector r is inside one of the atomic spheres and return to the second 
of our main equations, Eq. (|15ip . then we see that there are two possibilities: either the integration variable is on the 
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surface of the same sphere, or another one. For the first of these possibilities we have s, > 7"j and we can therefore 
write the Green function as 

G+(r,s)=Y(f l ) T p»(r l )(Mr + 1 ) T q (+)( Sj fY(s i ). 
The surface integral may be then evaluated as follows: 

/ [G+ (r, s) V S V; + (s) - (s) V S G+ (r, s)] • n t da, 



Y(f i ) T Pl (r l ) (M-;; 



Or, 



q l i+J ( Sl ) T Y(s l )Y(s l ) T R l '( Sl )C i 
da, 



-q| +) ( Sl ) T Y(s l )Y(s l ) T R l ( S ,)C l 
= b t 2 Y(r t ) T Pl (n) (M^f [qW(6 l ) T R i '(6 l ) - q i (+) '(6 l ) T R l (^) 
For the second possibility, with j ^ i, we have Tj > Sj, and the Green function is given by 

G+(r,s) = Y(r i ) T q; +3 (r i ) (MJ+ 1 ) P^) T Y(%) , 
which becomes, noting that rj < Rij and using the re-expansion formula (|139p . 

G+(r,s) = Y(f l ) T p l (r l ) (M^) T K « (M^ 1 ) Pi (^) T Y(%) . 

We then have 

/ [G+(r,s)V s ^+(s) - ^+(s)V s G+(r,s)] ■ n,- do-,,- 

JdTj 

= bfVptf Pi (n) (M^) T K « (M7+ 1 ) 'p,«/ ); ;'RM/ ); i p/:/,, i' R<;/,j 

If we define the Wronskian of two matrix functions a(r) and b(r), say, by 

W[a, b](r) = a(r)b'(r) - a'(r)b(r) , 

then our results may be written 
= 

X+(r) + Y(f i ) T p i (r i ) (M^f J &* 2 W[q( +)T , R l ](6,)C + £ 6/ K« 



(7 



i)W[ P /,R^)C7 



Let us now define 

= bi 2 (Mrl)\N[ Pi T ,R i ](b i )C i 

which we may presumably invert to recover 

C i = 6r 2 W[p/,R i ](6 i )" 1 M i+ B i 
Define also, for each atomic site, the matrices 



T! 



" 1 =W[q < (+)T > R i ](6 i )W[p i T J R i ](fc i )- 1 



(157) 



(158) 



(159) 
(160) 
(161) 

(162) 



We now note that the inhomogeneous term may be written, using Eqs. (73]) and (fTTHf (noting that the re-expansion 
we wish to effect is from the origin of coordinates to center j), as 



X +(r) = (4jt/A) Y(?. i ) T p. i (r. i ) (M^) D(*;Ri)£Y(k) . 



(163) 
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We may now write our results 



= Y(f,) T Pl (r l ) (M^f { (4w/k) D(fc; R t ) £Y(k) + T* +J2 ^ & \ ; (164) 



for which it is necessary that 



T a 1 B i + Yl K%j& = -(47rA)D(A:;Ri)^Y(k). (165) 



There will be one of these equations for each i, and together they constitute the multiple-scattering equations for 
the case of an incoming plane wave of direction k. Clearly they form a linear system which may be solved for the 
unknown vectors B 4 by inversion of the secular matrix whose elements are defined by 

C=T^,^ + (l-%)^- (166) 

More generally, following Ref. 0], we replace the exciting amplitude 47r«'Yj,(k) by Sll" and find a solution vector 
B l L (L") for each partial- wave channel L". The results may be then applied to more general situations. 

With the solution vector B l we may find the wavefunction inside the atomic spheres by using Eq. (|161[) to determine 
the coefficients C, then inserting them into Eq. (|152p . For the interstitial-region wavefunction we must return to 
perform the surface integrals in Eq. (|150p . In this case, for r in the interstitial region and s on the surface of atomic 
sphere i, we have r\ > Sj, so that 

G+(r,s) = Y(f l ) T q l (+) (r l ) (M^) p 4 ( Sl ) T Y(s,) . 

We obtain then 

N 

V>+(r) - x +(r)+^& 4 2 Y(r,) T q l (+) (r l )(M- 1 ) [p i (6 i ) T R i '(6 i ) - p/(6 i ) T R i (6 i ) 

z=l 
JV 

= X + (r)+^Y(r l ) T q| +) (r,)BV (167) 

i-l 

This completes the solution of the multiple-scattering problem. 



1. Symmetry of the secular matrix 



We see from Eq. (| 142[) that the second term in the secular matrix (|166p is symmetric. The question then arises as to 
whether the first term T l a 1 has any symmetry. There appears to be none manifest in its definition (|162[) . However, 
the author has treated (in Ref. [HI]) a similar problem for the bound-state version of the method of Natoli, Benfatto, 

and Doniach In fact, a proof of the symmetry of T* 1 may also be developed in our present case. 
To see this we first define the matrices d and t by 



d = W[p T ,R](&)- 



[Vi- 



and 



(d 



-l\T 



T„ 



(d- 1 ) 



suppressing for the moment the atomic index i. We may recover our original matrix via 

T.- 1 =d T td . 



It is clear from the definitions that the symmetry of t is equivalent to that of T a 
Expanding t we have (suppressing dependence on the radius b) 



W[p T , R] T ( 

{R' T [p(m 



1+ x ) T W[q 



(+)T ,Rl 



;Tq (+)J 



pWfq 



R > — f 



(168) 
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where we have defined 



f = R 



p(M^) T q«- 



R + R J 



p'(M; 1 ) T q(+) J 



R' 



(169) 



It may be seen now that the two expressions collected between braces in Eq. (|168[) are each individually symmetric; the 
first by virtue of Eq. ([26]) and the second from Eq. (f59|) . So the symmetry of f is equivalent to that of t. Substituting 
from Eq. (|27|) in the second term of Eq. (|169|) leads to 



R' 1 



p(m;T q W 



R+ R T 



(M;!)p T 



R' 



1 

62 



R T R' . 



The term in braces here is now manifestly symmetric and we are left with the term in R T R'. However this last is also 
symmetric for similar reasons to the comparable term in p of Eq. (|24[) . i.e., since R is the matrix solution, regular at 
the origin, of the coupled radial Schrodinger equations for a real potential. The proof may be seen in Ref. [l|. Thus 
our result is shown. 



2. Bound States 



As mentioned above our new MSW method is readily adapted to bound states. In their exposition of a non-MT 
MSW method, Natoli, Benfatto, and Doniach [5[ do not treat this case explicitly (this is done in Ref. [15(] ) , but observe 
that it is achieved by dropping the inhomogeneous term in the Lippmann-Schwinger equation; making the analytic 
continuation k — > ik, E — > — E where necessary; using the appropriate re-expansion formulas for modified spherical 
Bessel, Neumann, and Hankel functions; and imposing decaying-wave boundary conditions in the asymptotic region. 
In our case we expect that similar comments should apply. 

Although again we shall not develop the bound-state case explicitly, it is necessary to address briefly several issues 
which arise in our present approach. Firstly, we note that dropping the inhomogeneous term in the Lippmann- 
Schwinger equation (|146[) leads to 

V(r) = J G+(r, S )V A (s)^(s) d 3 s (170) 

as the main equation to be satisfied by a bound state tp(r) at some energy E < 0. It is not difficult to show directly 
that such a ip(r) will satisfy Schrodinger's equation for the full potential V(r), even though the two-stage picture of 
the distorted-wave formalism is lost. 

A more substantial problem arises when we come to derive the re-expansion formulas. In particular, since there is 
no inhomogeneous term in Eq. (|170[) . the argument leading to the re-expansion formula (| 1 1 8[) for the general regular 
solutions is no longer available to us. In principle, following the comments above, we may analytically continue 
Eq. (|118[) to obtain a similar equation: 

Y&fpifa) (MT'f = Yfcfp^-) (MJ^ D^( K ; Ry) ; (171) 

where k is such that k > and n 2 = —E : and D(~)(«;Ry) is obtained from D(/c;Ry) by the substitution k — > ire, 
therefore having components 

D { ~1 (k; R) = 4vr ^ i'-'i,,, (kR)Y l „ (R)7(L, L' , L") . (172) 

L" 

The Pi(fi) and qj(rj), etc., are solution matrices at our negative energy E, and the qj(r"i) are asymptotic to the 
decaying-wave free solutions q/(rj), given by 

q/(r)=i£k+ (/cr), (173) 

rather than those of Eq. (|4l))) . since this is the analytic continuation of the free positive-energy solutions (|28p. Such a 
procedure, however, should really be more rigorously mathematically justified. 

An alternative derivation may be developed via an argument more in line with our earlier approach. We first note 
that, taking p/ from Eq. (|4"4")l and q/ from Eq. (|173p . we find the inverse transpose of the constant matrix M/ to be 



(174) 
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The re-expansion formulas for the free solutions at negative energies are well known (see Eqs. 11.25 and 11.23 of 
Ref. 0] for example). They may be also obtained by the substitution k — > in in our earlier Eqs. (|9H)) and ([57| . It is 
not difficult to check that they can be written in our present notation as 

Y(r l ) T p/(r l ) [uf) T = Y(f,) T p / (r J ) (m?) T D<->(k;R«) (175) 

and 

Y(r i ) T q/(r i ) = Y(rj ) T qj (r 3 - ) D^ - ^ (k; Ry ) for r 3 - > R l3 . (176) 

If we now consider a situation such as that of Fig. [TJ with no restriction on r< other than that we maintain the order 
relations Sj > Rij, Sj > rj, and Sj > rj, then we may use the general Green- function formula to write, similarly to 
Eq. (TTT91) . 

Y^ftfcHMr 1 )^^)^^ (177) 

If, for the moment, we assume that the potential is identically zero beyond some large radius R, and if we take the 
sphere of possible vectors s, large enough to be entirely outside this radius, then c\j(sj) and qi(s^) may be taken to be 
identically q/(sj) and q/(si) respectively in this region. We may therefore use the transpose of Eq. (|1T6[) to substitute 
on the lhs of Eq. (|177p , rearranging then to obtain 

Y(r i ) T p i (r i ) (M-'f D^(k; Ryf - Yfofpjfo) (Mjfj q/ ( Sj ) T Y(s,) = , 

and we may conclude that the expression in square brackets here is identically zero. Now by an argument similar to 
that leading to Eq. ([95]) we may show that 



Dt-^R)- 1 = D(-)(k;-R) . (178) 

Furthermore it is easily seen from the explicit formula (|172|) that 

D(->(k;-R) = D(-)(«;R)* = D^^Rf (179) 

(showing incidentally that )(k; R) is hermitian), and we may therefore deduce that the re-expansion formula (|17ip 
for the general regular solutions at negative energies is indeed correct. Presumably we may use this result to obtain, 
in a similar way to the derivation of Eq. (|120[) , the corresponding general far-region irregular-solution re-expansion 
formula 

Y(s l ) T q,:(s t ) = Y(B,-) r q J -(«i)D(-)(/s;R ii ) for Sj > R l3 . (180) 

Our derivation here is less general than we should like because of the assumption that the potential is zero beyond some 
large radius. If, however, we note the fact that the regular solutions Pi(r^) are determined by the choice of boundary 
conditions at = 0, and that the D'~)(k;R) are completely given quantities, then the effect of the potential cutoff 
will only be felt by the q^(r^) and consequently the M^. Thus, since we may keep rj and rj fixed in Eq. (|17ip . only 
Mi and Mj will vary as a function of the potential cutoff, and we may plausibly invoke a limiting argument to justify 
the formula as the cutoff tends to infinitely large distance, i.e., no cutoff. 

A further issue arises when we consider the potential V/(r). It is possible (indeed likely for large molecules) that 
Vj(r) itself may support bound states. This would appear to cause problems when our energy passes through that of 
such a state. In the bound-state multiple-scattering case the secular equation is constructed in a similar way to that 
for the continuum case (|165p . except that there is no inhomogeneous part. The condition for a bound state is that 
the secular matrix becomes singular, allowing a non-zero solution vector from which the corresponding wavefunction 
may be constructed. As shown in Ref. [l| a bound state of Vj(r) leads to the constant matrices Mj, involved in the 
construction of the secular matrix, becoming themselves singular, and would give rise to problems when the inverses of 
these matrices are required. Conceivably it may be possible to factor out such a singularity, or deal with the problem 
in another way. This issue remains to be investigated in more detail. 

A systematic way of bypassing this problem is suggested by a comment of Natoli, Benfatto, and Doniach in the 
context of their non-MT MSW method Q. In particular, we combine their subtraction of a constant value from 
the interstitial-region potential (useful in their case to obtain a convergent Born expansion for the interstitial-region 
scattering T matrix) with the use of an outer sphere (as we describe in the next section) to obtain a Vj (r) which does 
not bind. 
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A final issue is that of normalization of the wavefunction. Since the wavefunction is constructed from the coefficients 
B l in the solution vector of the homogeneous bound-state secular equation, the normalization is not known. One must 
integrate the square of the resulting function over all space to determine the multiplying coefficient which leads to 
unit norm. In practice one may evaluate this numerically (although there are analytic methods that may be adapted 
to our present approach [4{) and we shall leave further consideration of this issue outside the scope of the present 
work. 



3. Outer Sphere 

As in previous versions of the MSW method it is sometimes useful to add to the partition of space the outer sphere 
to, containing all the atomic spheres, that we have defined in an earlier section. In the region to outside the outer 
sphere the potential is expanded as a multipole expansion relative to the sphere center Ro and the coupled radial 
Schrodinger equations are solved, similarly to the atomic centers, leading, for r 6 to, to 

^(r)=Y(f ) T R°(r )C , (181) 

for some unknown coefficients C , with the difference being that the R°(to) are integrated inwards from suitable 
boundary conditions at infinity. 

The use of an outer sphere leads to a modified secular equation. We shall not reproduce the algebra here since it is 
a relatively straightforward adaptation of the method of Ref . (BJ , using the formulas and approach we have developed 
above. 

The utility of the outer sphere is that it enables one to impose boundary conditions directly on the solution of the 
Schrodinger equation, important in the case of, for example, long-range Coulomb potentials. It is also useful to reduce 
the size of the interstitial region, which is important in the MT version of the MSW method, and in the non-MT 
version of Ref. [H, it helps to reduce the effect of the interstitial T matrix. For our present approach maximizing 
the size of the atomic spheres (and minimizing the size of the outer sphere) would tend to produce smaller potential 
gradients at the sphere surfaces and hence a more gently varying distorting potential Vj (if the atomic parts of Vj 
were chosen to be smooth continuations of the true interstitial potential) . 

As mentioned in the previous section, the use of an outer sphere is a way of avoiding the possibility that Vj has 
bound states of its own. This may be done by choosing the extra-molecular part of Vj to be asymptotic to some 
constant potential sufficently less than zero that no bound states arise in the resulting new Vi. Since the old Vj 
must be bounded from below, one obvious choice would be the minimum of this potential, although it may not be, 
in general, necessary to set the new asymptotic value as low as this. For particular cases some experimentation may 
be required. The smooth continuation of the new Vi into the extra-molecular region may be then set identical to 
the chosen asymptotic value beyond some large radius, so that the resulting q, ; would be identical to the diagonal 
modified spherical Hankel function forms for all larger values of the radius. 

One consequence of this new effective zero for the Vi is that between this value and the true energy zero of V we 
must use the standing- wave Green function related to Vj, in a manner similar to the original version 0] of the MSW 
method. This is, however, expected to be a straightforward modification of the usual positive-energy formulas (|22[) . 
with the c\i(ji) now asymptotic to n(k ri) and k calculated relative to the new asymptote of Vi- 



III. NUMERICAL IMPLEMENTATION 



The main computational task for the solution of the multicenter scattering problem is the setting up of the secular 
matrix which is then inverted to solve the secular equation. The essential ingredients for this are the radial solution 
matrices p iy q\ + \ and R l for each atomic center. For these we require the multipole components at each center of Vi 
from r,; = to oo, and of V from r*j = to 6j. These, however, need only be calculated once. The radial solutions 
must be calculated for each energy required. With these solutions we can evaluate the constant Wronskian matrices 
Mi + at n — hi, which are then inverted; and the matrices for which the surface integrals may also be evaluated 
at the same radius. From these requirements we see that and R* are needed over the range ri = to bi, and 
from ri = bi to oo. (Infinity here may be taken to be a suitably large radius.) 

The solution of coupled radial Schrodinger equations for various types of potential has received considerable attention 



over many years and has a sizeable literature (see for example Ref. 2J]), and, for our immediate purposes, we may 
consider currently available numerical methods adequate. Although it is known that there are situations where 
problems may arise, for the types of potential we consider they appear not to occur. Our own experience with the 
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matrix Numerov method, for example, in the form given in Ref. [5[, has thrown up no obvious difficulties; and, for 
calculations where an independent check is available, gives quite acceptable results. 

The key to an accurate implementation of our present approach is the choice of distorting potential Vi(r). It is 
important that this potential be well behaved, avoiding discontinuities and wide variations. Since it is constrained to 
coincide with the true interstitial potential in that region, the only freedom we have is in the choice of values inside 
the atomic volumes. Here we should replace the true singular atomic potential with a smooth continuation of the 
interstitial potential which is finite everywhere. This may be done in a variety of ways. One possibility, which has 
been used in a slightly different context by the author, is to replace, for rj < 6$, each multipole component of the true 
potential V by a polynomial, in r,-, of minimal degree, which matches the true multipole component at r% = and 
(perhaps) some of its derivatives. 

IV. DISCUSSION 

Our new approach to non-MT MSW theory results in formulas similar in overall form to those of some previous 
versions, and indeed should be thought of as a direct development of them. Although the derivation of the new 
secular equation in particular owes much to the approach of Ref. [f| , the final form of the result is more reminiscent 
of the original MT-based version of Slater and Johnson [2] . This may be understood if one notes that our distorting 
potential corresponds to the constant potential of their method. Thus where one sees free-electron waves (expressed 
in terms of spherical Bessel, Neumann, and Hankcl functions) in their formulas, one sees the distorted waves Pi and 

in ours. (In fact, the MT version of the method should result from our present algebra when Vj is replaced by the 
constant average interstitital potential.) We treat directly the multiple scattering in terms of waves which are already 
being scattered in the distorting potential. By contrast, and very loosely speaking, the method of Ref. [j| adds as an 
extra step (via the interstitial T matrix) the effect of the interstitial region potential on the waves scattering from the 
atomic centers. 

To enable development of the non-MT MSW theory in the direction we have taken it requires, however, elaboration of 
the formal machinery that makes up the bulk of this present work. It seems likely that some of the results derived here 
are new and should prove useful beyond our immediate subject. Firstly, the explicit general Green- function expansions 
of Ref. Q are necessary to implement the distorted-wave Lippmann-Schwinger formalism, leading, in particular, to 
an explicit form (|73|) for the distorted wave itself. Further direct development of the material of Ref. [l[ was required, 
leading to some useful formulas, notably the important relation (|59[) between regular and irregular matrix solutions 
of the coupled radial Schrodinger equations (|18p . We then saw how one can use these results to develop a complete 
(if skeletal) solution to the single-scattering problem for a general noncentral potential, whose most important result 
is the formula (|79p for the matrix A in the expression for the scattering amplitude. Indeed, having been presented 
with such a solution one might reasonably ask why it could not be applied directly to the scattering problem of the 
multicenter potential V. As is well known, the singularities of V have only slowly convergent spherical- harmonic 
expansions relative to a single center. To have an accurate treatment of scattering in such a potential requires an 
approach which is capable of combining individual accurate solutions in the regions of the singularities. This, of course, 
is the essential utility of the multiple-scattering approach to the solution of Schrodinger's equation for multicenter 
potentials. 

As mentioned in Ref. [l| certain mathematical issues relating to the derivation and manipulation of our Green- 
function expansions have been glossed over, and thus must be counted among the assumptions of our approach. The 
use of vectors and matrices of infinite dimension and questions of convergence should be considered in more depth. 
Such issues are frequently bypassed in the literature, it being noted that for practical purposes such vectors and 
matrices are truncated at finite dimension in real calculations. Of more physical interest is the asymptotic behavior 
of the general irregular radial solutions q;. The derivation of the Green-function expansions in Ref. 11 assumes that 
it is possible to choose these solutions asymptotic to a diagonal form. Indeed we assume here that this diagonal 
form is the free solution appropriate to outgoing waves of Eq. (|28|l . In the absence of finite charge moments the 
potential will typically decay exponentially and this would almost certainly be the case. It may easily be shown that 
the potential multipole expansion is dominated by the monopole component at infinity, even in the presence of a 
finite dipole or higher moment. Thus the potential matrix of Eq. (|19[) will be diagonally dominant. The nature of 
the limiting behavior necessary may in fact be seen from the expression (|7ip in the derivation of the explicit form for 
the scattered wave in a given potential, and in the algebra [at Eq. (|87p ] which verifies the defining equation for the 
general Green function. 

It is possible to bypass the issue of the asymptotic behavior of the irregular solutions by using the outer-sphere 
option. Here one may set the distorting potential to a constant value beyond some finite radius. In this region, as 
we mentioned above, the general irregular solutions may be chosen identical to the free ones, although once inside 
the finite radius they will, in general, deviate from them. This can be the case even if the true potential in the 



30 



extra-molecular region contains charge monopole, dipole or higher components. 

The second important set of results which enables our version of the non-MT MSW theory is that related to the 
re-expansion of general regular and irregular solutions of the coupled radial Schrodinger equations. This is central 
to our approach and we have devoted a significant fraction of this present work to a careful presentation of the 
derivations, returning more than once to particular points, and including alternative derivations of some formulas. 
As mentioned in the section on bound states, and illustrated in Eq. (171j) . the particular form of the matrix D(fc; Ry) 
associated with translations is dependent on the asymptotic form of q. It is interesting that this matrix, which has 
been encountered before in MSW theory, has a more general function than previously appreciated. 

V. CONCLUSIONS 

We have presented here a new approach to multiple-scattering theory for general potentials, based on the distorted- 
wave (two-potential) Lippmann-Schwinger formalism. (We propose to refer to it therefore as the "DW-MS" method.) 
Our method makes use of the familiar partition of space into non-overlapping atomic spheres and a remaining inter- 
stitial region; splits the multicenter molecular potential into a well behaved background distorting potential and a 
singular atomic part; and uses some newly developed Green- function expressions and re-expansion formulas to derive 
the new secular equations. 

The central improvement offered by this approach is to avoid any volume integrals, which will have significant 
benefits for the scaling behavior of calculations, the most complicated parts of which are now the surface integrals 
associated with the matrices K y . In fact, the possibility arises, perhaps by taking the limit of the integrals as 
Ti — > 0, that even these may be reduced to one-dimensional integrations. Calculation of the effects of the varying 
interstitial potential appears implicitly via the irregular solutions q[ + ^ of the coupled radial Schrodinger equations for 
the distorting potential Vj (r) . Another attractive feature of the approach is that the above calculational tasks require 
only numerical methods that are already well tried and tested in related approaches. 

It may be noted that, in a formal sense at least, the separation between geometry and scattering has been recovered. 
In reality of course this is not quite true since there is a connection between the distorting potential and the scattering 
centers. However, it may be possible to exploit this in an approximate way to reduce the computational burden in 
the context of geometry searches and total energy minimization. 
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FIG. 1: Diagram of vectors and sphere surfaces in the derivation of the re-expansion of the irregular radial solutions in the far 
region. 




FIG. 2: Diagram of vectors and sphere surfaces in the derivation of the re-expansion of the irregular radial solutions in the 
near region. 



